Direct observation of the degree of quantum correlations using photon-number 

resolving detectors 
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Optical parametric down-conversion is a common source for the generation of non-classical cor- 
related photonic states. Using a parametric down-conversion source and photon-number resolving 
detectors, we measure the two-mode photon-number distribution of up to 10 photons. By changing 
the heralded collection efficiency, we control the level of correlations between the two modes. Clear 
evidence for photon-number correlations are presented despite detector imperfections such as low 
detection efficiency and other distorting effects. Two criteria, derived directly from the raw data, are 
shown to be good measures for the degree of correlation. Additionally, using a fitting technique, we 
find a connection between the measured photon-number distribution and the degree of correlation 
of the reconstructed original two-mode state. These observations are only possible as a result of the 
detection of high photon number events. 

PACS numbers: 42.50.Ar, 42.50.Dv, 42.65.Lm 



Non-classical states of light are an essential resource 
for novel protocols in quantum information and quan- 
tum metrology 0, Q • The most common tool for produc- 
ing such states is the nonlinear process of optical para- 
metric down-conversion (PDC). In this process, a parent 
pump photon is split in a nonlinear material into two 
daughter down-converted photons, while conserving en- 
ergy and momentum. As the down-converted photons 
originate from a single quantum system, they possess 
correlations in many degrees of freedom, such as their 
polarization, frequency and momentum. Because for any 
photon emitted into one optical mode there is a sister 
photon emitted into the other optical mode, there are 
also photon-number correlations between the two modes. 
These correlations have been used to produce heralded 
Fock states [3mJ and enhance the precision of optical 
measurements [J, [2, 0, 0] ■ 

The down-converted photons are distributed over a 
range of spatial and spectral modes. However, most 
experiments require that the photons occupy a single 
mode. A specific mode is then post-selected by spa- 
tial and spectral filtering of the photons. The collected 
modes must be carefully matched to obtain a high-quality 
produced state jjij. The collection of matching modes 
would result in a joint photon-number distribution of a 
non-classical highly-correlated state. On the other hand, 
collecting two unrelated modes would result in a clas- 
sical joint photon-number distribution which is a prod- 
uct of the two individual states. Recent developments in 
photon-number resolving detectors allow the direct mea- 
surement of the joint photon-number distribution and 
photon-number correlations between two down-converted 
modes. However, imperfections in the detection pro- 
cess, such as collection losses and false detections, alter 
the measured photon statistics and reduce their correla- 
tions. Non-classical correlations of multimode distribu- 
tions from PDC have been demonstrated indirectly using 



reconstruction techniques lllLll2l and directly with a sys- 
tem of relatively low loss [ljj LLJ| 

In this Letter, we use a novel photon-number resolving 
detection scheme in order to measure the joint photon- 
number distribution of a collinear type-II PDC process. 
These distributions were measured up to the 10 photon 
terms. The two down-converted photons, which have or- 
thogonal polarizations, are each collected from a single 
spatial and spectral mode. The degree of correlation was 
controlled by varying the amount of overlap between the 
two collected modes. We directly observe the transition 
between a separable product state and a highly correlated 
state, despite the presence of low detection efficiency and 
other distorting effects. We introduce measures for the 
degree of correlations between the two modes and relate 
these measures to the degree of non-classicality of the 
collected state. 

The PDC states are generated by a type-II collinear 
(5 — BaB204 nonlinear crystal. The crystal is pumped by 
amplified and frequency doubled ThSapphire laser pulses 
at a repetition rate of 250 kHz. The orthogonally polar- 
ized down-converted photons at 780 nm are split using a 
polarizing beam-splitter and coupled into separate silicon 
photomultiplier (SiPM) photon-number resolving detec- 
tors [Hj (Hamamatsu Photonics, S10362-11-100U). Be- 
fore coupling to the detectors, the down-converted pho- 
tons are spatially filtered using single-mode fibers and 
spectrally filtered using 3 nm bandpass filters, to ensure 
the collection of a single spatio-temporal mode [17| ■ The 
degree of correlation g of the measured photon state is 
determined by the amount of spectral and spatial over- 
lap between the two collected modes. The amount of 
overlap is tuned by translating the optical fibers in or- 
der to collect different spatial modes and by tilting the 
bandpass filter, in order to shift its spectral band. The 
amount of overlap is evaluated using the heralded effi- 
ciency 7, defined as the ratio between the coincidence 
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Figure 1. (color online) (a) Photon number resolution of a 
single detector, demonstrated by the good peak separation of 
the pulse intensity histogram, (b) The joint photon-number 
distribution of two polarization modes (H and V) of a type-II 
collinear PDC process for a maximally correlated state with 
7 = 0.18. 



and the single count rates in the limit when the average 
number of photons approaches zero. The parameter 7 is 
linearly proportional to the degree of correlation between 
the modes g and to the overall photon detection efficiency 
r\. 7 was evaluated using standard photon-number non- 
discriminating detectors (Perkin Elmer SPCM-AQ4C). 

The SiPMs are composed of a two-dimensional array of 
avalanche photo diodes (APD) operating in Geiger mode. 
These detectors are operable at room temperature and 
require a relatively low operating voltage. They gen- 
erate a current which is proportional to the number of 
detected photons. The output signals are sampled simul- 
taneously within a 1 ns sampling window, and analyzed 
in real-time using programable electronics. A computer 
receives the results and continuously displays the joint 
photon-number distribution between the two polariza- 
tion modes. The dark count and afterpulsing rates (l6| . 
which generally limit the number of resolvable photons, 
are minimized by synchronizing the sampling time of the 
SiPM analog output with the arrival time of the photons. 
Furthermore, the detectors are moderately cooled using 
a thermoelectric cooler to ~ — 10°C. The bias voltage 
is adjusted accordingly, so that the detection efficiency 
is not affected. As a result, we obtain a good photon- 
number resolution with an error of less than 1%. Figure 
[T][a) shows an example of the histogram of the output 
intensities, in which up to 14 photons can be resolved. 

A measurement of the joint photon-number distribu- 
tion between the horizontal and vertical polarization 
modes with maximal spatio-temporal overlap is shown 
in Fig. mb). We measured the full joint probability ma- 
trix up to 10-photon terms. Ideally, the photon-number 
distribution of the correlated state would be composed 
of only diagonal elements which correspond to the same 
number of photons in both modes. However, the distri- 
bution exhibits a large number of non-zero probability 
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Figure 2. (color online) Joint distribution measurements be- 
tween (a) highly correlated modes (7 = 0.18, blue circles) and 
(b) non-correlated modes (7 = 0.06, green up triangles). The 
corresponding products of the two individual distributions are 
shown as black down triangles. Poissonian errors are assumed 
and presented when larger than the symbol size. 



values for events which contain different photon num- 
bers. Similar probability values were measured for events 
which involve the same number of photons, indicated by 
the different color groups in Fig. [IJb). Note however 
that the probability of events which contain no photons 
in either one of the modes are slightly higher than the 
remaining probabilities within the same group. 

The distortion in the photon-number distribution is 
a result of the imperfect detection process in SiPM de- 
tectors [17j ■ Not every photon impinging on the de- 
tector creates a signal due to imperfect overall photon 
detection efficiency. Additionally, false signals can be 
generated by thermally-excited discharges (dark counts) . 
Furthermore, when a detection element is triggered, it 
might trigger additional neighboring elements due to op- 
tical crosstalk [16], in which a spurious photon gener- 
ated during a discharge in one APD element propagates 
and is detected by another element. These inherent ef- 
fects alter the photon number distributions, reducing the 
photon-number correlations even between highly corre- 
lated modes. 

In order to study the photon-number correlations in 
our measurements, we compare the measured correlated 
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distribution with its corresponding non-correlated prod- 
uct distribution, obtained by multiplying the photon- 
number probabilities P(n^) and P(n v ) of the individual 
polarization modes, (n v ) is the number of photons in 
the horizontal (vertical) polarization mode. This result 
is shown in Fig. HJa) . The distribution is displayed such 
that each curve connects events which contain the same 
total number of photons S = ri/j + n v , and the event 
counts are presented as a function of the photon number 
difference D = — n v . This representation is similar to 
that of the joint distribution in Fig. [2(b), if each color 
group were to be joined with a solid line. The corre- 
lated distribution shows a clear distinction from its cor- 
responding product result. On the other hand, a similar 
measurement between two relatively uncorrelated modes, 
results in a distribution which is almost identical to its 
corresponding product result (see Fig. [U(b)). 

We have recorded a series of joint probability distri- 
butions for different values of the heralded efficiency 7. 
By maintaining a fixed value for the overall detection 
efficiency, the heralded efficiency can be used MS H di- 
rect measure of the mode overlap. In order to quantify 
the deviation of the joint distribution from an uncor- 
related product state, we define the ratio between the 
joint photon-number probabilities and the corresponding 
product of the probabilities of their individual polariza- 
tion modes 

R(n h ,n v ) = P(n h ,n v )/P(n h ) ■ P(n v ). (1) 

As can be seen from Fig. [Ha), the values of R are 
fairly uniform for all probability values except for the ex- 
tremes, when there are zero photons in one of the modes. 
The average values of R obtained for all probabilities 
R(nh 7^ 0, n v 7^ 0) for the series of distributions are pre- 
sented in Fig. |31[a) . These values approach R sa 1 for the 
least correlated states and increase linearly up to R ss 2 
for the highly correlated state. Numerical calculations of 
the dependence of R on the degree of non-classicality in 
the original state confirmed the linear relation between 
the two parameters. The exact value of R depends on 
the specific values of the average number of photons in 
the original distribution, the detection probability, the 
dark count rate and the crosstalk probability. Thus, 
even though it may seem that the original correlations 
are completely washed out, the ratio between the joint 
and product probabilities can be used as a direct measure 
for the degree of correlation between the two modes. 

A mathematical tool which provides a more universal 
quantitative measure for the similarity between a given 
two-mode distribution and its closest product state is the 
singular value decomposition [l8[ ■ The decomposition of 
an n x n probability matrix results in n singular values 
Si, ordered as si > S2 > ■ ■ ■ > s n > and normalized 
to J2i s i = 1- The decomposition of a separable product 
state results in a single non-zero value, si 7^ 0, whereas a 
maximally correlated state results in n values of Si = 
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Figure 3. (color online) Quantitative measures for the de- 
gree of correlation of the joint probability distributions, (a) 
The ratio R of the joint probabilities to their respective prod- 
uct values, (b) The Euclidean distance ||AM|| between the 
measured distribution and its closest product state, (c) The 
reconstructed degree of non-classicality g. 



1/y/n. The normalized Euclidean distance between a 
measured matrix M and the closest product state M pr d 
can be expressed using the singular values of M and is 
given by ||AM|| = \\M - M prd \\ = + s 2 3 + ■ ■ ■ + s 2 n . 

The Euclidean distances ||AM|| for several heralded 
efficiency values are shown in Fig. ^h). The errors in 
||AAf|| are estimated using a bootstrapping procedure 
and assuming Poissonian noise. The distance exhibits a 
clear increase by an order of magnitude as the heralded 
efficiency is varied between its minimal and maximal val- 
ues. Thus, the Euclidean distance 1 1 AM| | is a measure for 
the degree of correlation of the two-mode state. It is ap- 
plied directly to the raw data, and can detect correlations 
despite large imperfections in the detection apparatus. 

In the experiment, the heralded efficiency 7 was first 
optimized such that the largest amount of coincident 
events were observed. This condition is fulfilled when 
the two collected polarization modes maximally over- 
lap, both spatially and spectrally. Then, by misaligning 
one of the collected modes, the overlap between the two 
modes was reduced, as well as the photon-number cor- 
relations. Thus, the joint probability of the two down- 
converted modes is a linear combination of the probabil- 
ity Ppdc of a correlated PDC distribution and that of an 
uncorrelated product distribution P pr d 

P(n hl n v ,g) = g-P p dc(n h ,n v ) + (l-g)-P prd (n h ,n v ) , (2) 

where the degree of correlation g is the amount of over- 
lap between the two collected modes. For g — 1, the 
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Figure 4. (color online) Measurements of distributions with 
fixed product statistics and different heralded efficiencies of 
7 = 0.17 (clue circles), 7 = 0.11 (green diamonds), 7 = 0.02 
(red triangles) and the calculated product state (black tri- 
angles) . The solid lines are fits to Eq. [2] with the following 
parameters: average number of photons (n) = 4.1±0.1, detec- 
tion efficiencies: r) A = 0.012 ±0.005 and -q B = 0.010 ± 0.0005, 
crosstalk probabilities: e A = 0.12±0.01 and e B = 0.11±0.01, 
average dark count rate: Aa = 0.11 ± 0.01 and As = 
0.14 ±0.01. 



distribution is that of a collinear type-II PDC 

'0 n h ^n v 



P p dc(nh,n v ) 
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n>h = n v = n . 



(3) 



where the parameter (n) is the average number of pho- 
tons in each mode. For g = 0, the distribution is that of 
a product state of two thermal modes 
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By modelling the effects of loss, dark counts and 
crosstalk in the SiPM detectors, we reconstructed the 
original two-mode photon-number distribution using a 
fitting method, similar to that presented in Ref. [l7|- The 
measured distribution P m is related to the original dis- 
tribution P as 



Pm 



(5) 



where P m and P are vector representations of the mea- 
sured and the original two-mode probability matrices, re- 
spectively. The matrices M; oss , Mdk, and M C ( represent 
respectively the effects of loss, dark counts and crosstalk 
according to the detector model of Ref. [I?} • The fitting 
procedure is performed in two stages. First, we perform 
a least squares fit to the computed product state and ob- 
tain the average number of photons per mode (n), the 
overall detection probability 77, the average number of 
dark counts, and the crosstalk probability for each of the 
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Figure 5. (color online) The reconstructed photon-number 
probability distributions for the data presented in Fig. H ( a 
7 = 0.17, g = 0.47 (b) 7 = 0.11, g = 0.23 (c) 7 = 0.02, 
g = 0.06 and (d) the product state. 



two modes. Then, the measured data is fitted to Eq. [SJ 
with the degree of correlation g as the free parameter. 

Some example results of the fitting process are shown 
in Fig. 2] Three distributions taken for different val- 
ues of 7 are presented. The photon-number probabilities 
in both polarization modes were kept constant for all 
measurements, thus maintaining the same product state 
(black lines) . Transition from a highly correlated state to 
a product state is observed as the value of the heralded 
efficiency is decreased. The solid lines in Fig.@]are fits to 
Eqs.[2] All fits result in similar values for the dark counts, 
detection efficiency, crosstalk probability, and (n) up to 
the experimental error. The fits clearly differ in their 
values for the degree of correlation g. 

Figure [5] shows the reconstructed distributions from 
the data of Fig. 0J For the highest correlated state we 
observe strong photon-number correlations, which grad- 
ually disappear as the value of the heralded efficiency 7 
is decreased. Even for the lowest measured value of 7, 
photon- number correlations are still evident. 

The degree of correlation g for additional values of 
7 are shown in Fig. [3tc). As expected, the degree of 
correlation g depends linearly on the heralded efficiency. 
The parameter g is also the degree of non-classicality 
of the state. Applying the non-classicality criterion for 
the photon statistics of two-mode radiation of Lee [ijj on 
Eq. [2] shows that a necessary (but not sufficient) condi- 
tion for non-classicality is that g > 0. We have tested all 
of the reconstructed distributions against Lee's criterion 
and found that all states with g > satisfy it. 

In conclusion, we have measured the two-mode photon- 
number distribution of a collinear type-II PDC process 
for different degrees of correlation. Clear evidence for 
photon-number correlations are presented despite the 
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low detection efficiency, the dark counts and the opti- 
cal crosstalk effects of SiPM number-resolving detectors, 
which highly distort the number-correlations. These ob- 
servations are only possible as a result of the detection 
of high photon number events. Both the ratio between 
the measured and the product matrix probabilities, and 
the singular value decomposition of the measured prob- 
ability matrices, are shown to be good measures for the 
degree of correlation. These two criteria are derived di- 
rectly from the raw data. Additionally, using a fitting 
technique, we found a connection between the measured 
photon-number statistics and the degree of correlation of 
the reconstructed original two-mode state. 

The authors thank O. Gat for fruitful discussions. 



[1] P. R. Tapster, S. F. Seward, and J. G. Rarity, Phys. Rev. 

A 44, 3266 (1991). 
[2] J. Dowling, Contemp. Phys. 49, 125 (2008). 
[3] C. K. Hong and L. Mandel, Phys. Rev. Lett. 56, 58 

(1986). 

[4] E. Waks, E. Diamanti, and Y. Yamamoto, New J. Phys. 
8, 4 (2006). 

[5] P. J. Mosley, J. S. Lundeen, B. J. Smith, P. Wasylczyk, 
A. B. U'Ren, C. Silberhorn, and I. A. Walmsley, Phys. 
Rev. Lett. 100, 133601 (2008). 

[6] J. Rarity, P. Tapster, and E. Jakeman, Opt. Commun. 



62, 201 (1987). 
[7] H. S. Eisenberg, J. F. Hodelin, G. Khoury, and 

D. Bouwmeester, Phys. Rev. Lett. 94, 090502 (2005). 
[8] S. F. Seward, P. R. Tapster, J. G. Walker, and J. G. Rar- 
ity, Quantum Optics: Journal of the European Optical 
Society Part B 3, 201 (1991). 

[9] W. Mauerer and C. Silberhorn, Phys. Rev. A 75, 050305 
(2007). 

[10] W. Wasilewski, C. Radzewicz, R. Frankowski, and K. Ba- 

naszek, Phys. Rev. A 78, 033831 (2008). 
[11] E. Waks, E. Diamanti, B. C. Sanders, S. D. Bartlett, and 

Y. Yamamoto, Phys. Rev. Lett. 92, 113602 (2004). 
[12] M. Avenhaus, H. B. Coldenstrodt-Ronge, K. Laiho, 

W. Mauerer, I. A. Walmsley, and C. Silberhorn, Phys. 

Rev. Lett. 101, 053601 (2008). 
[13] E. Waks, B. C. Sanders, E. Diamanti, and Y. Yamamoto, 

Phys. Rev. A 73, 033814 (2006). 
[14] A. Allevi, A. Andreoni, F. A. Beduini, M. Bondani, M. G. 

Genoni, S. Olivares, and M. G. A. Paris, EPL 92, 20007 

(2010). 

[15] G. Bondarenko, B. Dolgoshein, V. Golovin, A. Ilyin, 
R. Klanner, and E. Popova, Nucl. Phys. B. 61, 347 
(1998). 

[16] P. Buzhan, B. Dolgoshein, L. Filatov, A. Ilyin, V. Kaplin, 
A. Karakash, S. Klemin, R. Mirzoyan, A. Otte, 

E. Popova, et al., Nucl. Instrum. Meth. A 567, 78 (2006). 
[17] L. Dovrat, M. Bakstein, D. Istrati, A. Shaham, and 

H.S. Eisenberg, Opt. Express 20, 2266 (2012). 
[18] E. Henry and J. Hofrichte, Methods Enzymol. 210, 129 
(1992). 

[19] C. T. Lee, Phys. Rev. A 42, 1608 (1990). 



